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ABSTRACT 

Type la supernovae (SNe la) are currently the best probes of the dark energy in the universe. To constrain 
the nature of dark energy in a model-independent manner, we allow the density of dark energy, px(z), to be 
an arbitrary function of redshift. Using simulated data from a space-based supernova pencil beam survey, we 
find that by optimizing the number of parameters used to parametrize the dimensionless dark energy density, 
f(z) = Px(z)/px(z = 0), we can obtain an unbiased estimate of both f(z) and the fractional matter density of the 
universe fl m (assuming a flat universe and that the weak energy condition is satisfied). A plausible supernova 
pencil beam survey (with a square degree field of view and for an observational duration of one year) can yield 
about 2000 SNe la with < z < 2 (Wang 2000). Such a survey in space would yield SN peak luminosities with a 
combined intrinsic and observational dispersion of cr(m„„) = 0.16 mag. We find that for such an idealized survey, 
fi m can be measured to 10% accuracy, and the dark energy density can be estimated to ~ 20% to z ~ 1.5, and ~ 
20-40% to z ~ 2, depending on the time dependence of the true dark energy density. Dark energy densities which 
vary more slowly can be more accurately measured. For the anticipated SNAP mission, Q m can be measured 
to 14% accuracy, and the dark energy density can be estimated to ~ 20% to z ~ 1.2, Our results suggest that 
SNAP may gain much sensitivity to the time-dependence of the dark energy density and 0,„ by devoting more 
observational time to the central pencil beam fields to obtain more SNe la at z > 1 .2. 

We use both maximum likelihood analysis and Monte Carlo (when appropriate) to determine the errors of 
estimated parameters. We find that Monte Carlo analysis gives a more accurate estimate of the dark energy 
density than the maximum likelihood analysis. 
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1. INTRODUCTION 

Supernova data suggest that most of the energy in the uni- 
verse is unknown to us (|Garnavich et al. 1998aj R iess et 
al. 1998, |Perlmutter et al. 1999[ ). Much theoretical effort 
has been devoted to exploring t he possible candidates f or the 
dark energy (for example, see Peebles & Ratra 1988L F rie- 
man et al. 1 995, Caldwell, Dave, & Steinhardt 19981 S ahni 
& Wang 2000), and investigating the constraints on the nature 
of the dark energ y that can be derived from observation al data 
(for example, see White 19981 iGarnavich et al. 1998bj S tein- 
hard t, Wang, & Zlatev 1999, pfstathiou 1999L Huterer& Turner 
2000, podariu & Ratra 20001, |Waga & Frieman 2000i Ng & 
Wiltshire 2 001, Podariu, Nugent, & Ratra 200 U W eller & Al- 
brecht 2001). 

The most straightforward and promising probe of dark en- 
ergy is the distance-redshift relation derived from observa- 
tions of cosmological standard candles. At present, type la 
supernovae (SNe la) are our best c andidates for cosmologi- 
cal s tandard candles (Phillips 1993,_Riess, Pres s, & Kirshner 
1995, |Phillips et al. 19991 pranch et al. 200l[). Most work- 



ers have concentrated on constraining the equation of state of 
the dark energy wx from SN data. However, without making 
specific assumptions about wx (for example, assuming it to be 
a constant), it is extremely hard to constrain wx using SN dat a 
( |Maor, Brustein, & Steinhardt 20001 parger & Marfatia 2001] ). 
Recently, Wang and Garnavich (2001) [WG01] showed that it 
is easier to extract constraints on the dark energy density px(z), 
instead of wx(z), from data. This is because there are multi- 
ple integrals relating Wx(z) to the luminosity distance of SNe 
la, diXi), which results in a "smearing" that obscures the dif- 

'wang@mail.nhn.ou.edu 



fere nce between different wx(z) (M aor, Brustein, & Steinhardt 
2000). On the other hand, px(z) is related to the time deriva- 
tive of the comoving distance to SNe la, r'(z); hence it is less 
affected by the smearing effect. The advantage of measuring 
px(z) over that of measuring wx(z) is confirmed by Tegmark 
(2001). 

WG01 gave a proof of principle of an adaptive iteration 
technique for extracting the dimensionless dark energy density 
f(z) = px(z)/px(z = 0) as an arbitrary function from future SN 
data. They found that feasible future SN data will allow us to 
clearly differentiate dark energy with density that changes with 
time from a cosmological constant A; however, estimates of ft,,, 
and f(z) tend to be significantly biased for dark energy densities 
that vary substantially with time. 

In this paper, we apply a significantly improved and opti- 
mized version of the adaptive iteration technique for obtaining 
an unbiased estimate of the dark energy density from SN la 
data. 

2. DARK ENERGY DENSITY 

The total energy density of the universe is 



p(z) = pl(l+z? + p° k {l+zj 1 +Pxf(z) 

= p° c [O m (l + z) 3 + n k (l+ zf + n x f(z)] 



(1) 



where the superscript "0" indicates present values, f(z = 0) = 1, 
and = 1 — il m — ilx = —kjH§. If the unknown energy is 
due to a cosmological constant A, f(z) = 1. Clearly, f(z) = 
Px(z)/px(z = 0) is a very good probe of the nature of the un- 
known energy. Following WG01, we impose the weak energy 
condition (the energy density is nonnegative for any observer 
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[ |Wald 198J ), which implies that f'(z) > (see WG01). 
The equation of state of the dark energy is 



Px(z) 3 f(z) 



1. 



(2) 



A constant wx(z) corresponds to f(z) oc (1 +z) a , where a is a 
constant. The values a = 0, a = 3, and a = 4 correspond to a 
cosmological constant, matter, and radiation respectively. We 
parametrize f(z) as 



f(z)= , 

K Zi~Zi-\ 

Zo — 0, Zn — Zmax 



fi-i + f ' 1 ) /i, Zi-1 < Z < Zi, 



^Zi~ Zi-i 
fo = l 



(3) 



where /) (i= 1,2, ...,n— 1) are independent variables to be esti- 
mated from data. We let /„ be either an independent variable, 
or linearly extrapolated from /„_i and /„_2, whichever gives the 
smaller \ 2 P er degree of freedom for the same n. 

The measured distance modulus for a SN la is /jfP = 
fiy + e (/) , where /i® is the theoretical prediction fjp = 

5 log (rfi(z/)/Mpc) +25, and e (/) is the uncertainty in the mea- 
surement, including observational errors and intrinsic scatter in 
the SN la absolute magnitudes. Denoting all the parameters to 
be fitted as s, we can estimate s using a modified \ 2 statistic, 
which results from integrating the probability density function 
for parameters s, p(s) oc exp(~x 2 /2), over the Hubble constant 
H Q . We write (see WG01) 
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where 



A.* 



Ci = 



*(/) (0 



C 2 = 



(4) 



(5) 



(6) 



j u; = /i i ,(ft = r) = 42.384-51ogr+51og[ffoKl+z)]- (7) 

We take ft* = 0.65. Our results are independent of the choice of 
h*. 

The current cosmic microwave background (CMB) 
anisotropy m easurements seem to indicate that we liv e in a 
flat universe ( |de Bernardis et al. 2000] ftalbi et al. 2000| ). Clus- 
t er abundances strongly suggest a low matter density universe 
( iBahcall, Lubin, & Dorman 1995j|Carlberg et al. 1996|, B ahcall 



& Fan 1998). fl m = 0.3 and tt A = 0.7 is the best fit model to 
current observational data. We will use fl m = 0.3 and fix = 0.7 
for our simulated data in the rest of this paper. 

In order to compare with the results of WG01, we consider 
the same two hypothetical dark energy models, given by 



Mz) = 



,1.5c 



/,(z) = exp[0.9(l-O], 



w 9 (z) = -l+0.5z 

w k (z) = Q3(l+z)e- z -l; (8) 



fq(z) and fk{z) represent quintessence (dw g /dz > 0) and k- 
essence (dwt/d z < 0) models respectively (C aldwell, Dave, & 
Stein hardt 1998, A rmendariz-Picon, Mukhanov, & Steinhardt 
2000) and they satisfy the weak energy condition f'(z) > 0. 



3. RESULTS 

A feasible SN pencil beam survey (with a square degree field 
of view and for an effective observational p eriod of one year) 
can yield almost 2000 SNe la out to z = 2 ( |Wang 2000"| ). Let 
us combine the data from the SN pencil beam survey with SN 
data at smaller redshifts, so that there are a minimum of 50 SNe 
la per 0. 1 redshift interval at any redshift. This yields a total of 
1966 SNe la for quintessence and 1898 SNe la for the k-essence 
model, up to z = 2 and for Q m = 0.3, fix = 0.7. We assume in- 
trinsic and observational dispersions which are Gaussian with 
zero mean and a variance of 0. 16 magnitudes (for a space-based 
survey). 

Fig. 1 shows the likelihood function £(0,,,) oc <rx n ™( n <*)/ 2 for 
the quintessence model f q {z) and the k-essence model fkiz). 
Note that Xmin^m) i s marginalized over n independent param- 
eters that parametrize the dimensionless dark energy density 
f(z), fi (i = 1,2, ...,ti) (see Eq.(^|)). For the quintessence model 
fq(z), we allow /„ to be an independent variable, and for the 
k-essence model fk(z), we linearly extrapolate /„_2 and f n -\ to 
obtain /„ (see text after Eq. [[?]]). The curves (peak location from 
left to right) correspond to n = 3,4,5, 10 and n = 4,5,6, 10 for 
f q (z) and f k (z) respectively. 
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FIG. 1. — The likelihood function £(Q m ) oc exp [-xL„(^m)/2] for 
the quintessence model f q (z) and k-essence model fkiz). 

For the quintessence model f q {z), the \ 2 P er degree of free- 
dom decreases as we decrease n from n = 10 to n = 4, with a 
modest shift in the most likely value of fl m (see Fig. la); this is 
as expected because our estimate of Cl m should be roughly in- 
dependent of the parametrization of f(z). There is a substantial 
shift in the most likely value of Vl m as we change n from n = 4 
to n = 3. Although the \ 2 P er degree of freedom is smaller for 
n = 3, the optimal choice is n = 4 for this model, since a greater 
degree of degeneracy between f2 m and /, (i = 1,2, ...,«) sets in 
for n = 3, which renders the estimated value of Q m significantly 
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biased. For real data, and the true value of f2 m unknown, this 
sudden increase in degeneracy can be inferred via the substan- 
tial shift in the estimated value of fl m . 

For the k-essence model fkiz), the % 2 per degree of freedom 
decreases as we decrease n from n = 10 to n = 4, with the largest 
shift in the most likely value of f2 m ocurring as we change n 
from n = 5 to n = 4 (see Fig. lb). Hence n = 5 is the optimal 
choice for this model. Note that the transition that marks the 
sudden increase in degeneracy (at n = 4) is less dramatic than in 
the case of the quintessence model f q (z). 

Fig. 2 shows the dimensionless dark energy density f q (z) and 
fit(z) estimated with n = 4 and n = 5 respectively (see Fig. 1). 
The dashed line and circles represent the estimate derived for a n 
ideal space-based supernova pencil beam survey (Wang 2000). 
The dotted lines and triangles represent the estimate derived for 
the SNAP 2 mission. The errors are (16%, 84%) intervals of the 
cumulative probability distribution of the estimated parameter 
for 10 3 Monte Carlo samples. 
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FIG. 2. — The dimensionless dark energy density f q (z) and ft(z) es- 
timated with n = 4 and n = 5 respectively (see Fig.l). The dashed line 
and circles represent the estim ate derived f or an ideal space-based 
supernova pencil beam survey ( |Wang 2000| ). The dotted lines and 
triangles represent the estimate derived for the SNAP mission. 

Note that Fig. 1 and Fig. 2 differ significantly from Fig. 2 
of WG01, where the estimate of both ft m and /, (i =1,2, ...,«) 
become significantly more biased as n is changed from n = 10 
to n = 6. This work is greatly improved over that of WG01 as 
follows: (1) Maximum likelihood analysis is used here to de- 
termine the optimal choice of n in parametrizing f(z), which 
was not explored by WG01; (2) The end point of f(z), /„, is 
allowed to differ from f„-\ here, while /„ = /„_i was imposed in 



WG01; (3) The grid sizes of Afi m = 0.001, A/, = 0.05 are used 
here, compared with the grid sizes of Ail n , = 0.02, A/i = 0.1 of 
WG01; (4) The errors of the estimated parameters are derived 
via Monte Carlo from 10 3 instead of 10 2 random samples as in 
WG01. 

We have used a combination of maximum likelihood analy- 
sis and Monte Carlo technique (when appropriate) in this paper. 
Fig. 3 and Fig. 4 show the probability distribution functions of 
fl m and (i = 1,2,3,4) for the quintessence model / 9 (z). The 
dotted curve is the likelihood function derived from marginal- 
ized x 2 , while the histogram is the probability distribution de- 
rived from Monte Carlo of 10 3 random samples. The arrows in 
each figure indicates the true value of the parameter. Clearly, 
Monte Carlo analysis gives more accurate estimate of the dark 
energy density f(z) than the maximum likelihood analysis for 
large f(z). 
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FIG. 3. — The probability distribution functions of Q,„ for the 
quintessence model f q (z)- The dotted curve is the likelihood function 
derived from marginalized x 2 , while the histogram is the probability 
distribution derived from Monte Carlo of 10 3 random samples. 

4. DISCUSSION 

Due to a substantial improvement in our adaptive iteration 
technique (see WG01), we are able to explore much greater 
ranges of possibilities in the parametrization of the dimension- 
less dark energy density f(z) = px(z)/px(z = 0). As a result, we 
have developed a method for finding unbiased estimates of O m , 
together with f(z) = px(z)/px(z = 0) (as an arbitrary function 
parametrized by its values at n equally spaced redshifts). The 
optimal choice of n corresponds to the smallest x 2 per degree of 
freedom without significant shift in the estimated bestfit value 
of il m (see Fig.l). This leads to the unbiased estimate of fi m , 
which is crucial in deriving an unbiased estimate of the dimen- 
sionless dark energy density f(z). 

Our approach is more robust and model-independent than 
constraining the dark energy equation of state wx(z) (as a linear 
func tion of z) assuming that fl,„ is known (H uterer & Turner 
2000, |Weller & Albrecht 2001[ ), and complementary to the lat- 
ter in probing the nature of the dark energy. 

We find t hat for an ideal supernova pencil beam survey 
(Wang 200C) from space, f2 m can be measured to 10% accu- 
racy, and f(z) = px(z)/px(z = 0) can be estimated to ~ 20% 
to z ~ 1.5, and 20-40% to z ~ 2, depending on the time de- 
pendence of the true px(z) (see Fig. 2). Dark energy densities 
which vary more slowly can be more accurately measured, as 
might have been expected, since it is easier in general to con- 
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Fig. 4. — The probability distribution functions of /, (i = 1 , 2. 3, 4) for the quintessence model f q (z) with n = 4. The line types are the same as in 
Fig. 3. 



strain slowly varying functions compared with rapidly varying 
functions. 

For the anticipated SNAP mission, il,„ can be measured to 
14% accuracy, and f(z) = px(z)/px(z = 0) can be estimated to 
~ 20% to z ~ 1.2 (se e Fig. 2). Co mpared with the idealized SN 
pencil beam survey ( Wang 2000[ ), the SNAP strategy is to ob- 
tain a larger number of SNe la (2408 versus 2000 for a 1° x 1° 
pencil beam) at z < 1 .2 by devoting a significant fraction of the 
observational time on flanking 1 x 1 fields that surround the 
two central 1 x 1 ° pencil beam fields, at the price of a sharply 
decreased number of SNe la at z > 1.2; this leads to large er- 



rors in estimated /,■ for z > 1 .2 (see Fig. 2). Our results suggest 
that SNAP may gain much sensitivity to the time-dependence 
of px(z) and fi m by devoting more observational time to the 
central pencil beam fields to obtain more SNe la at z > 1 .2. 

Future supernova pencil beam surveys hold great promise for 
constraining the nature of dark energy in a model-independant 
manner. 

We thank Eddie Baron and David Branch for helpful com- 
ments. YW acknowledges support from NSF CAREER grant 
AST-0094335. 
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